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In this work, the cosmological imphcations of brane world scenario are investigated when the 
gravitational couphng G and the cosmological term A are not constant but rather there are time 
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^ H and A ~ H^. The behavior of scale factors and different kinematical parameters are 
investigated for different possible scenarios where the bulk cosmological constant A5 can be zero, 
positive or negative. 
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I. INTRODUCTION 

There are two fundamental physical parameters namely the gravitational coupling G 
and the cosmological term A, in the Einstein's theory of general relativity. Usually these 
parameters are assumed to be constants. But recent experimental and (or) observational 
evidences suggest that it is not unlikely to have small variations [1] of these fundamental 
parameters. However, alternative ideas about the variability of these parameters have been 
started long ago. The idea of variable G was first introduced by Dirac [2] in cosmological 
consideration and then subsequently by Sciama [3], Jordan [4] and others. Later, Brans 
and Dicke [5] proposed an extension (or modification) of Einstein's theory of gravity by 
introducing a scalar field (p. In this theory, the Newton's gravitational constant G is a 
variable and is related to the scalar field cj) ^ G~^. 

The cosmological term A on the otherhand, was originally suggested by Einstein himself 
to obtain static solution of his field equations. But after his realization that the universe 
is expanding, he discarded the cosmological term. Afterwards, the cosmological term has 
been introduced and subsequently rejected several times for various reasons. But recent 
measurements of the CMB anisotropy and the observations from type la supernovae demand 
a significant and positive cosmological constant [6, 7]. Also observations of gravitational 
lensing indicate the presence of a non-zero A. However, to resolve the cosmological constant 
problem (i.e., its observed value is about 120 orders of magnitude below the value for the 
vacuum energy density predicted by quantum field theory) a phenomenological solution 
namely a dynamical A{t) has been suggested by several authors [8], arguing that A relaxes 
its present estimate due to the expansion of the universe. Also due to these recent obser- 
vational predictions the vacuum energy density pA {— A/SttG) has played an important 
role and the cosmological constant problem has been shifted to the coincidence problem 
namely why p and pA happen to be of the same order of magnitude precisely at this very 
moment. To resolve this naturalness several quintessence models [9] have been proposed 
where the cosmological constant becomes a time dependent quantity. Recently Bonanno et 
al [10] (see also Shapiro et al [11]) have formulated a general frame work for cosmologies in 
which both Newton's constant and the cosmological constant are time dependent and have 
shown that p and pA are approximately equal in the late universe. They have assumed 
that there exists an infrared attractive fixed point for the renormalization group flow of the 
(dimensionless) Newton's constant and cosmological constant respectively. Vishwakarma 
[12] has proposed a variable A proportional to [H — Hubble parameter) and has showed 
a good agreement with recent supernovae observations. Also this choice provides expected 
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large age of the universe. 

In this paper, we analyse a 5D branc world model with a time dependent brane tension 
A and constant 5D cosmological constant A5. We shall show that this will lead to a 4D 
effective theory with time dependent G and A. Although the fundamental theory is not yet 
known, the time dependent of the self-energy of the brane A can be the result of some higher 
dimensional effect, such as a particular bulk-brane interaction, or the response of the brane 
to the change in the inter-brane distance in a 2-brane world model. Thus here we attempt 
to make a "phenomenological" study of the implications of such a time-dependence of A 
on observable 4-dimensional parameters, such as G and A, and its consistency with current 
observational constraints. 



II. BASIC EQUATIONS IN BRANE COSMOLOGY 

The idea of brane world scenario was proposed by Randall and Sandrum [13]. They have 
shown that in a five dimensional space-time (called bulk) it is possible to confine the matter 
fields in a four dimensional hypersurface (called 3-brane). The effective action in the bulk 
is [14] 

A = J d^X^^^ (2^^^"^^)"^ / ^'^^^V^ (^K^ C-^atte^ (1) 

Here any quantity in the bulk is marked with subscript '5' with ^7^^) as the 5D 

gravitational coupling constant. The co-ordinates on the bulk arc denoted by x". a = 
0, 1, 2, 3, 4 while those on the brane as , = 0, 1, 2, 3 . So x = a;^ = the four dimensional 
hypersurface, is defined as the brane world. The parameter A stands for the brane tension 
which is assiimcd to bo positive to recover the conventional gravity on the brane. The other 
parameter A5 is the negative vacuum energy and is the only source of gravitational field on 
the bulk while are the intrinsic curvature on either side of the brane (characterized by 
X > or < 0) . Now the effective Einstein equations on the bulk are [15, 16] 



<^a6 — "^S^ab + '^5^a6(brane) (2) 

where T'^b(brano) (with T'ff,(brane) ~ 0) total energy-momentum tensor (vac- 

uum-|-matter) on the brane i.e.. 



7^a6(bra„c) = ^{x) [Affafc + T^^\ = KxWa^r^u ■ (3) 

Here T^^^ is usual matter energy tensor and the Dirac delta function indicates the fact 
that matter is confined in the space-like hypersurface x = (the 3-brane) with induced 
metric gab- 

From the Israel's junction conditions 

Due to Z2-symmetry of the bulk, 

Ki^u\x>o = -Kij,^\x<o = ^nl (ji^u - ^9ij.i.T ) (4) 

Hence we obtain 

2 

'^5 
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So from (4) we have 



Ki_iu\x>0 — ~2'^5 



^^17^ - + A) 
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Then the effective Einstein equations on the brane are ([15]; see also [17]) (choosing 
87r= l,c= 1) 



G^^ = -Ag^, + K^T^™) + KjSi,. - E^^ (5) 

(Note that Greek indices refer to brane world while English small letters to bulk quantities). 

The above Einstein equation has two correction terms to the energy-momentum tensor 
namely the local correction term S^v with expression 

C _ ''" rp(m)rp(m) ^ rpa(m)rp(m) i "'^ „ /orpaf3{m,)rp{m) rp2(m,)\ /f>\ 

while the other correction term is the non-local effects from the free gravitational field in 
the bulk with expression 



Eab = CabcdU^n'^ (7) 

Here Cabcd is the usual Weyl tensor in the bulk, n°' the unit normal to the hypersurface 
X = and the four dimensional cosmological constant A has the expression 

A = ^(A5 + ^) (8) 

with 

n'' = ^ (9) 

as the four dimensional gravitational constant. Now using the conservation of the total 
matter energy-momentum tensor on the brane namely 

. = (10) 

into the modified Einstein equations on the brane we have a constraint on 5^,^ and Ei^^, as 

{El-4Sl),, = Q (11) 

Maartens has shown that E^i, can be decomposed as [16, 17] 

E^,u = — ^[U{Uf,u, + ^ V) + 2w(^(3.) + P^A (12) 

K A o 

4 

where the prefactor (i.e., ^) is introduced for dimensional reasons (note that A'^ ^ 
gives the general relativistic limit) and the energy flux and energy stresses P^i, has the 
following properties : 

(i) Q/u is a spatial vector (i.e., Q/nW'' = 0) 

(ii) Pjui/ is a spatial (i.e., P^„u" — 0), symmetric {P{fj.v) = Pfii^) and trace-free {Pj^ ~ 0) 
tensor. 
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The scalar term U is called dark energy density as it has the same form as the energy 

momentum tensor of a radiation perfect fluid. Also is the 4-vclocity on the branc with 
hfj_^, the projection tensor orthogonal to u'^ on the brane. Further, using the constraint 
equation (11) between and S^i,, it is possible to have evolution equations for U and 
but not for So the system of equations on the brane is not, in general closed. 

As we arc considering Bianchi models which have a simple transitive three dimensional 
group of isomctries G3 on the space-like hypersurfaces, so the cosmological models are 
spatially homogeneous with proper time as the only dynamical variables. Thus the metric 
on the brane can be written as [18] 

ds^ = -dt'^ + a^i{t)dx'dx' (13) 

Further, the physical parameters which are of observational interest in cosmology are the 
following: 



Expansion scalar: e = 3H = Hi = X;-=i ^ 

Shear scalar: 

1 „., 3 



Deceleration parameter: q = — ^^h^ ^ 

where the mean anisotropy parameter A has the expression 



and Hi and H are the directional Hubble parameters and mean Hubble parameters 
respectively. 

Lastly, for Bianchi I and V model we get = but P^^ remains unrestricted. As there 
is no way of fixing the dynamics of this tensor, we shall study the particular case P^i, = 0. 
Then from the constraint equation (11) the dark energy density C/ is a function of the proper 
time alone i.e., U = U{t) and the evolution equation for U takes the form 

U = - 4HU (15) 

which on integration gives 

U = UoV-i (16) 
where V = 11?=! '^i is the volume scale factor and Uq is an integration constant. 

III. BIANCHI I MODEL 

The line element of a Bianchi I space-times generalizes flat FRLW metric to the anisotropic 
case as 

rfs^ = -dt^ + al{t)dx'^ + al{t)dy'^ + al{t)dz^ (17) 
We have taken perfect fluid as the matter on the brane with expression 
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where the energy density p and thermodynamic pressure p satisfy the isothermal equation 
of state namely 

p = (7 - l)p, 1 < 7 < 2. (19) 

The combinations of the non-vanishing components of the Einstein field equations on the 
brane (i.e. (5)) are [18] 

3H + ±H! = A- (^^.V - ^^4p' + ^ (20) 

and 

It is to be noted that so far we have not considered any variation of the 4D physical 
parameters namely G and A. The variation of G and A shows a time variation of brane 
tension A (see equations (8) and (9)), keeping 5D quantities as invariant. Prom observational 
point of view it is generally assumed that the time variation of G can be expressed in terms 
of the Hubble parameter H as [19, 20] 

^=59H i.e., G = GoV^ 

where the dimensionless parameter g has present observational bound \g\ < 0.1 .and Go 
is constant of integration. 

As from equation (9) G oc A so the differential equation in A has the same form namely, 

J = 3gH i.e., A = Aq^s (22) 
where Aq is integration constant. Also the integration constants (Go and Aq) are connected 

4 

by the relation (using eq.(9)) Gq = ^Ao- 

As a consequence, the energy conservation equation (10) becomes 

p + 3jpH = -XogV^-^V (23) 

which on integration gives 

p = (24) 

with Ao = FqD and 1 + (3 = —g = 

Hence from observation /? is restricted to — 1 < /? < —0.966. Also from equation (9) the 
brane tension is proportional to G and we write 

X=^hzl^V-^^-^^ (25) 

Now using relations (8), (16), (22) and (23) in the field equation (21) we have 

Vdt^^^''~ 2 12(1 + (3)^^ ^ ^ 

So adding equation (26) for i = 1, 2, 3 we get 
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which on integration gives 



H = Hi + 



V 



(28) 



The integration constants Ki, i = 1,2,3 are restricted by the relation Ki = 0. 



Combining equation (26) with (27) we have the differential equation in V (after integrating 
once) as (assuming (1 + /?) 7^ 0) 



or in integral form 



dV 



(29) 



(30) 



where a = |k|A5, b = ^^f^) c = 9Uo and ((ii,fo) are integration constants. 

Also the scale factors can be obtained integrating once equation (28) and using equation 
(29) as 



fli = aioVs exp 



3 



dz 



, i = 1,2,3 



with Uio's as integration constants and z = V ^ . Further, if we assume the variation of the 
cosmological term (as mentioned earlier) as 



(31) 



(with ^, a function of time) then using (25) and (31) in equation (8) we get a differential 
equation in V as 
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y2 = ^^KjA^V^ + ^KtD^ 



3 



(7-1-/3)^ 
(l + /3)2 



,y-2/3 



Now comparing with equation (29) we have 



' (7-l-/3f 



+ cV^2(l+/3)+^^y2/3^0 



(32) 



which is true for all time and for all values of V. As powers of V in the 3rd and 4th term 

arc in no way compariable to the first two terms for any choice of /? so wc can choose the 
integration constant di = and also c — 0. Wc note that the constant 'c' is related to 
the non-local energy correction term (corresponds to an effective radiation) which is con- 
strained to be small enough at the time of nucleosynthesis and it should be negligible today . 

Moreover as a solution of (32) we have (in addition to c = = rfi) three possibilities 
namely 



{i) ^ = 3 and 7^ = (7 - 1 - pf i.e. 7 = 



(1 + /3) 



(33) 



(ii) a = and ^ 



3(7-1-/?)' 



(34) 



and 



(iii) ^ is a function of V. 
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It is to be noted that the first choice contradicts the observational bounds for /3 [19] (i.e., 

— 1 < /3 < —0.966), so only the second and the third possibilities will be discussed in the 
following for all possible values of A5 namely, (i) A5 > 0, (ii) A5 < and (iii) A5 = 0. 



Case I: A5 (i.e., a) > 

Here we have the solution of V from equation (30) as (the third possibility i.e., ^ is a 
function ofV, is the only choice) 



= ^1 sinh[(l + p)^ {t - to)] 
which simplifies to (choosing to = 0, so the big bang singularity occurs at t = 0) 



(35) 



with 



3(7-/3-1)2 ^ 4A, 



r 



+ 



2H^ 



The expressions for different physical parameters are 



(36) 



H=l = ^^ Kr, coth[(l +/3)^2A5 Kg 

P = cosech[yi|r «5(1 + 0) t] 

q = 3{l + p) sech'^[\f^ ^5(1 + (3) t] - 1 
A 



(7-/3-i)76AF cosech[^/M^ «5(1 + /3) t] 



a = 



-./EM 



V- 



^ = 3(7-/?-!)- g2 + [1 _ 



G = ^1(2^^ /^ cosech[./i|^ K5(l + /3) 



(37) 



From this solution we observe that H decreases sharply with time for small t and 

approaches a constant value >i5\J~^- This assures the positivity of p, A and G. Asymptot- 
ically, for large t, the universe expands exponentially with positive acceleration q^i —\ and 
A and G become constant (see figs. 1 and 2) while universe isotropizes due to exponential 
fall off of the anisotropic scalar (see fig. 3). Also the interdependence of fi^ and /3 over q 
has been shown in the 3D graph (fig. 4). 

Further, the expressions for density parameter and the age of the universe are 



2(7-/3-l)(l + g) 
372 
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the acceleration or deceleration of the universe at different cosmic time for A5 > 0. 



si' 




Fig.3 Fig.4 



Fig. 3 shows the gradual isotropization of the universe with time. Fig. 4 gives the variation of q 
with the variation of Q.p and /3 for A5 / 0. 



tanli 



2 + 3/3 - g 



3(1 + /3) 



1/2 



So we can write 



^ ^ 2(1 + g) 

i.e., (3 can be estimated from the observed values. 
Again from the relations (37) we have 



Also from (30) using (37) we write 



7v/3(l + /3) 



H 



3 2, 3g^(2 + 3/3-g) 
2 ^ 1 + /3 



(38) 



(39) 



(40) 



(41) 



Hence from the above two relations we can evaluate the five dimensional quantities K5 and 
A5. Thus the solution contains no free parameter. 



We note that a non-vanishing cosmological constant in the bulk induces a natural time 
scale in AD as 
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which in terms of observed quantities has the expression 

^'=h[^- 2(7-/3-l)(l + /3)j (^^^ 

If we consider a model universe with flp « 0.3 and that of H 0.7 x 10"^'^ y~^, then the 
characteristic scale is » 5.2 x 10^° yrS i.e., 52 billion years, which is much greater than the 
age of the universe (due to data from WMAP mission) [21]. 

Lastly, from the observational bound for (3 the deceleration parameter is restricted by 
_ ]^ <; g <; 2{-y^034) ^ ^' whlch shows that for flp < ^^'^^^^^^^ , there is always an 
accelerating universe while for fip > ^^^^a'^^'' both acceleration and deceleration is possible 
(see fig. 4), but for ftp > only we have decelerating universe. 

Thus in the present matter dominated era (7 = 1) to get accelerated universe ftp < 0.644 
which supports the observed value of Op. Therefore it is possible to have at present an 
accelerating universe in this model. 



Ccise II: A5 (i.e., o) = 

Here the relation (34) should be satisfied and V has the simple solution 

1 

V = \Vb{p + l)t\ , (/? + 1 ^ 0) (43) 

with 

3(7-1-/?)^ 

(Note that 0+1^0; otherwise there is no variation of G ) 

Here also the integration constant to is taken to be zero for simplicity. The scale factors 
will have an explicit forms as 



tti = ttioV'^ exp 



, i = l,2,3 



3/3 

The deceleration parameter will have constant value and is related to /? by the relation 

g = 2 + 3/3 

The brane tension, AD effective cosmological term A and the gravitational coupling G 
has the time variation as 

Ki7(/3 + l)i' 7' 7 ^ ' 



with 
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q 




Fig.5 Fig.6 

Figs. 5 and 6 present the variation of q due to the variation of f2p and 7 for '+' and ' — ' sign in 
equation (46). 



as Hubble parameter. 

The anisotropy scalar a decreases with time as t^T+^ and vanishes at infinity. Now for 
energy density namely p — , we have the following relation among density parameter 

rip, /3 and 7 

Gp 2(7-/3-l)(/3 + l) 

2 = "p = -5 (45) 



which gives /9 as a function of 7 and Qp 

7 



P + l = j[l± v/1 - 2f^p] (46) 

Thus for real f3 we must have ilp < 1/2. This upper limit of Qp is independent of any 
specific values of 7 and q. Also eliminating between (45) and the expression for q we 
have 



The variation of g over f2p and 7 has been shown in figures 5 and 6 for positive and 
negative signs in equation (46). In the first case there is always deceleration while in the 
second case q varies between —1 and 2. 

Now the relationship between the density parameters are 

= "a = \, (48) 



3i/2 4(1 + g)2 

Further we note that the above solution can be recovered from the solution for A5 > 
if we proceed to the limit A5 0. Near the big bang epoch (i.e., for small t) the present 
solution has similar behaviour with that for positive A5. Lastly, to have an accelerated 
universe the density parameter should have further restriction on upper bound (depending 
on 7), considering sign in equation (46) as 
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We note that this upper bound for Op is always less that 0.5, hence it is possible to have 

a decelerating phase of the universe when Qp lies between the bound in equation (49) and 
0.5. On the other hand, for '+' sign in equation (46) the universe will accelerate if 7 < 2/3 
and then there is a lower bound for ftp i.e.. 



^^^^<n,<0.5 (50) 

To have an end of this section we remark that similar to the previous solution for A5 > 0, 
here also there is no free parameter remain, all of them can be estimated from observation 
(i.e.. Hp and H). In table I and II we have presented the numerical values of different 
parameters based on various values of the density parameter Qp for dust model (7 = 1) 
with H = 0.7 X 10~^° yr~^. Table I corresponds to '— ' sign in the expression for /? in 
equation (46) while table II for '+' sign. In table I, throughout there is acceleration except 
for flp — 0.5 as expected from the expression for q. In fact, for Hp > 0.44 only deceleration 
is possible. We note that here flp + Ha 7^ 1, there is a contribution of Hp2 in brane scenario. 
But this contribution is insignificant in the present epoch. Thus the data in this table are 
consistent with recent observation as there is deceleration followed by acceleration. On the 
otherhand, table II has the reverse picture. There is always deceleration as 7 = 1 and flp2 
has significant contribution in the matter content of the universe. Finally, we say that table 
II represents early phases of the universe when brane scenario has dominant contribution 
while table I is valid in the present epoch when brane effect is negligible. 



TABLE-I 

(A5 = 0. Present phase of the evolution, accelerated expansion) 



Hp 


P 


q 


Oa 


T/IO"' yr's 


0.5 


-0.5 


0.5 


0.25 


0.95238 


0.4 


-0.72361 


-0.17082 


0.52361 


1.72287 


0.3 


-().8i(i2;] 


-0.44808 


().()()()2;] 


2..')!)12() 


0.2 


-0.88730 


-0.66189 


0.7873 


4.22523 


0.1 


-0.94721 


-0.84164 


0.89721 


9.02108 


0.08 


-0.95826 


-0.87477 


0.91826 


11.4078 


0.04 


-0.97958 


-0.93875 


0.95958 


23.3234 
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TABLE-II 

(A5 = 0. Early phase of the evolution, decelerated expansion) 



ftp 


P 


q 




T/W yr's 


0.5 


-0.5 


0.5 


0.25 


0.95238 


0.4 


-0.27639 


1.17082 


0.07639 


0.65808 


0.3 


-0.18377 


1.44868 


0.03377 


0.58340 


0.2 


-0.11270 


1.66190 


0.01270 


0.53668 


0.1 


-0.05279 


1.84164 


0.00279 


0.50273 


0.08 


-0.04174 


1.87477 


0.00174 


0.49693 


0.04 


-0.02042 


1.93875 


0.00042 


0.48612 



Ccise III: A5 (i.e., a) < 

This case is important from the point of view of brane scenario as proposed by Randall 
and Sundrum [8] because according to them a single brane is embedded in a five dimensional 
anti-de-Sitter bulk (AdSs). 

For this choice the solution of V can be written as (with to = 0) 



= sin[(l+/3)v^t] 



(51) 



with the expressions for the physical parameters as 

= I = cot[(l + P)^/\a\ t] 

■s/elAsI (1+13) r,^ , „x /]— I 

P = ^ cosec[(l + P)y/\a\ t] 

q = 3(1+13) sec2[(l + /3)^t] - 1 

A 



(7-/3-1)^/6|a7| 



K57 



cosec[(l + /3)v|a| t] 



A 



3(7-/3-1)" 



(7-/?-!)" 



(52) 



Q = 4h-0-^) ^lJM cosec[(l + t] 

This solution represents a spatially flat but recoUapsing universe with recoUapsing time 
(-Tree) given by 



Sm[{l + /3)^/\^ Tree] =0 
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Though this solution is completely different from the solu tion with positive A5 but we can 

obtain this solution from (37) if we change y/A^ z-^/jAsj. As a consequence, \/a — > i-\/faf 
and sinh(-ya t) ^ i sm{^/\a\ t). Finally, the age of the universe can be written as 

1/2 



Tree ^ _i M - 2 - 3/3 

T = tan — — -— 

TT V 3(1 +/3) 



IV. BIANCHI V MODEL 

The Bianchi V space-time is an isotropic generalization of the open {k = —1) Robertson- 
Walker geometry and the metric-ansatz has the expression 



ds^ = -de + al{t)dx'' + al{t)e-^''dy^ + al{t)e-^'' dz'' 
As before the field equations are equation (20) and 



(53) 



^^.^p-t_^4p^ + ^,i = l,2,3 (54) 



2Hi - i?2 - i?3 = (55) 

Now compare to the Bianchi I model in the previous section, we have the same relation 
(28) among the Hubble parameters but the constants are given by 

Ki = 0, K2 = -Ks (56) 

due to the field equation (46). The differential equation for V in the present model is of 
the form (after integrating once) 



(57) 



Also as before considering the variation of A as A = (,H^, we have the same differential 
equation (34) as in Bianchi I model. Comparing (48) with (34) we note that Bianchi V 
model is possible only for /? = i (with c = = rfi) which contradicts the observational 
bounds on (3. The explicit solution is (choosing the integration constants appropriately) 



' 3^cosh-iyi^, if ^<o 

3^sinh-iyi^^, ff ^>0 



3^l0gy3/2, if i, = 4 



The explicit form for the physical parameters are : 



(58) 



H=^ tanh(|v^ t) 

p = D.f^ (cosh|VS i)-^ 



g = -1 - I cosech^di/a t) 



A=|D(7-i)y^(cosh|^r' 
A = Ci?2 



G oc (cosh|v^ i)~i 



> for < 



(59) 
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-1 



> for 



b-4 



>0 



(60) 




a 



A=|£»(7-i)J^ (sinhf^i) 



-1 



G oc (sinh|-ya t) 



-1 




g = -l 

A= §£'(7- I) e-i^* 



In this paper, we have studied the consequences of the variation of G and A in the form 
^ ~ 7? and A ~ on Bianchi I and V cosmological models based on the branc-world 
scenario. We have obtained solutions in Bianchi I model for positive, zero and negative 
value of the cosmological constant A5 on the bulk. Prom the solution with posit ive A5 
one can obtain solution for negative A5 by analytic continuation (v^As i-y/jAsJ) while 
solution for A5 = arc the limiting solution of (37) as A5 0. However, asymptotically, 
as <: — > 00 the nature of the solutions in three cases arc distinct. For A5 > 0, the universe 
expands exponentially with constant A as f — > 00, while for A5 = the universe expands in 
a power law fashion (for (3^—1) with a constant deceleration and all physical parameters 
become vanishingly small. But for A5 < 0, the solution is an expanding-collapsing model of 
the universe i.e., an oscillatory model with finite time period. 

For A5 > 0, and G fall off sharply from large positive value for small t to a constant 
value for large t, while A approaches to zero asymptotically. Also for A5 = 0, the behaviour 
of these parameters are similar but all of them vanish asymptotically. But we have peculiar 
behaviour of these parameters for negative A5 as there are points of discontinuities at finite 
t. Wc have also presented solutions for Bianchi V model with A5 > only for different 
signs (including zero) of the constant , but they arc not realistic as they corresponds to 
(3 = 0.5, which is outside the observational bounds on (3. Their mathematical form is very 
similar to the Bianchi I case, with one correction term, proportional to V^^^ added to the 
parametric time equation. Hence the general physical behaviour in both geometries has 
the similar qualitative features. Lastly, we have shown the behaviour of the shear scalar 
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(cr cx t~^+T) against t and /? in fig. 7 for standard cosmology Bianchi I model. In standard 
model a falls off as polynomial in 't' while in brane scenario it falls exponentially. In fact 
the figures 3 and 7 show the similarity of the anisotropy scalar in brane scenario and in 
standard cosmological model. 

Finally, we like to emphasize that the whole analysis of this paper is independent of any 
particular choice of the bulk parameters which is a good characteristic for the brane-world 
models. 
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